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Abstract 

This  paper  demonstrates  the  current  and 
future  potential  of  f inite-difference  methods  for 
solving  real  rotor  problems  which  now  rely  largely 
on  empiricism.  The  demonstration  consists  of  a 
simple  means  of  combining  existing  finite- 
difference,  integral,  and  comprehensive  loads 
codes  to  predict  real  transonic  rotor  flows. 

These  computations  are  performed  for  hover  and 
high-advance-ratio  flight.  Comparisons  are  made 
with  experimental  pressure  data.  .  ,  .,  > 
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*  xVc,  none  i  mens  ional  streamwise  coordinate 
(prime  denotes  a  physical  quantity) 

-  y’/R,  nondimensional  spanwise  coordinate 
(prime  denotes  a  physical  quantity) 

-  zvc«1/3,  scaled  nondimensional  normal 
coordinate  (prime  denotes  a  physical 
quantity) 

*  angle  of  attack,  deg 

»  partial  angle  of  attack,  deg 

-  specific  heat  ratio 

=  thickness  ratio 

-  collective  pitch  angle,  deg 

*  longitudinal  cyclic  input,  deg 

-  lateral  cyclic  input,  deg 

-  advance  ratio 

?  /•) 

-  j>’/r.Rc6  ,  a  scaled  velocity  perturbation 

potential 

*  normal  derivative  if  o 

*  second  normal  derivative  of  > 

*  azimuthal  angle,  deg 

=  angular  speed 

-  solidity 


Introduction 

The  design  of  new  rotor  systems  can  be  a 
-isky  undertaking.  Therefore,  most  of  these 
machines  are  derivatives  or  scaled  versions  of 
other  well-proved  rotor  systems.  Much  of  the  -isk 
of  helicopter  development  stems  from  an  inadequate 
understanding  of  and  ability  to  predict  the  per¬ 
formance-defining  physical  behavior  of  the  rotoi — 
especially  the  aerodynamic  behavior.  Most  of  the 
various  flow  phenomena  that  limit  or  define  the 
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capabilities  of  a  helicopter  (e.g.,  wake  forma¬ 
tion,  dynamic  stall,  acoustic  propagation,  tran¬ 
sonic  behavior,  blade/wake  interactions,  and  other 
unsteady  loading  and  drag  sources)  are  problems 
the  understanding  of  which  has  not  progressed  far 
beyond  the  empirical  level.  The  translation  of 
our  empirical  understanding  into  a  more  reliable 
predictive  ability  probably  awaits  the  development 
of  more  comprehensive  numerical  modeling  tech¬ 
niques  than  currently  exist.  At  present,  our  best 
computational  schemes  are  built  around  various 
integral  methods,  of  which  lifting-line  theory  is 
the  most  commonly  used.  These  methods  are,  in  the 
main,  limited  to  the  prediction  of  linear  compres¬ 
sible  flows.  The  use  of  tabular  airfoil  data 
extends  the  validity  of  these  methods,  but  only  to 
the  extent  that  two-dimensional,  steady  flow  is 
representative  of  the  rotor  flow. 

A  much  more  general  predictive  tool  is  to  be 
found  in  the  differential  methods  of  flow  model¬ 
ing.  This  class  of  schemes  (which  Includes  the 
f ini te-dif f erence  methods)  requires  a  high  invest¬ 
ment  t.n  computer  capability  but  has  no  intrinsic 
physical  modeling  limitations.  Finite-difference 
methods  for  the  solution  of  various  high-speed  and 
high-lift  problems  are  becoming  commonplace  in 
fixed-wing  aircraft  design.  However,  such  methods 
have  not  yet  been  employed  for  rotary-wing  uork 
because  of  the  greater  flow  complexity.  The 
intent  of  this  paper  is  to  demonstrate  a  means 
whereby  existing  codes  (both  integral  and  differ¬ 
ential)  may  be  combined  to  solve  realistic  rotor 
flow  problems,  which  have  heretofore  been  treated 
empirically. 


Choice  of  Flow  Models  for  Various  Flow  Problems 

The  helicopter's  various  flow  problems  can 
require  a  wide  number  of  models.  If  one  were  to 
compute  an  entire  vehicle  with  accurate  modeling 
of  all  -elevant  phenomena  nothing  less  than  a 
massive  ilavier-Stokes  (NS)  solver  would  be 
required.  But  far  more  modest  goals  have  con¬ 
founded  us  for  years.  The  oldest  of  these  is  the 
-otor  free-wake  problem.  All  integral  (vortex 
lattice,  etc.)  methods  have  failed  to  produce  a 
reliably  workable  tool — and  for  unknown  reasons. 

A  few  preliminary  Tlnlte-dlf ference  attacks1 '2  on 
the  problem  have  not  done  better,  but  this  work  is 
only  in  its  infancy.  The  best  current  approaches 
are  still  those  that  use  prescribed  wake  models. 
Retreating-blade  stall  is  another  problem  that  is 
best  handled  empirically.  The  problem  has  been 
treated  with  NS  codes,  but  these  are  not  yet 
satisfactory  because  of  Inadequate  turbulence 
model ing.  One  area  which  seems  no  longer  to  be 
beyond  our  grasp  is  the  high-speed,  advancing- 
blade  problem.  Separation  doe3  not  usually  occur 
here,  and  the  essential  pnysias  is  inviscld. 


Euler  codes  are  appropriate  for  these  flows,  but 
may  only  be  required  when  shocks  occur  that  gen¬ 
erate  appreciable  vorticity.  Flows  with  weaker 
shocks  are  much  more  efficiently  treated  by  poten¬ 
tial  methods.  Even  with  this  simplest  of  finite- 
difference  models,  it  is  not  currently  feasible  to 
treat  an  entire  rotor  flow,  and  the  problem 
becomes  one  of  efficiently  embedding  a  local 
finite-difference  solution  in  a  global  solution 
obtained  by  other  means  (probably  an  integral 
method).  Such  a  scheme  will  still  be  nonlinear, 
three  dimensional,  and  unsteady  and  could  greatly 
reduce  or  eliminate  the  need  for  airfoil  tables. 


Potential  Equations  of  Motion  and  Boundary 
Conditions 


The  starting  point  for  the  various  potential 
formulations  is  the  mass  conservation  and 
Bernouli's  equation,  together  with  the  assumption 
of  irrotationality .  This  system  of  equations  will 
conserve  mass,  energy,  and  entropy  when  solved  in 
a  divergence  form^  that  maintains  the  separate 
Identities  of  density  and  potential.  However,  a 
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number  of  significant  nonconservative  methods 
solve  only  for  potential.  Such  schemes  would  only 
be  conservative  if  combined  with  a  3hock-f ltting 
scheme.  (Note  that  conservation  implies  the  abil¬ 
ity  to  accurately  predict  shock  motion  and  the 
resulting  load  variations.)  The  simplest  poten¬ 
tial  equations  involve  various  small-perturbation 
approximations.  Thi3  class  of  equations  all  have 
the  form 
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where  F  is  a  nonlinear  r^ux  term  that  is  a  func¬ 
tion  of  *x-  This  equation  is  derived  by  means  of 
a  small  cross-flow  assumption  and  only  one  nonlin¬ 
ear  term  (that  contained  in  F)  is  retained.  A 
fortunate  consequence  of  the  small-perturbation 
approximation  is  that  the  equation  remains  conser¬ 
vative  in  3pite  of  it3  not  explicitly  involving 
density. 


The  boundary  conditions  for  Eq.  (1)  are  shown 
in  Fig.  t.  The  computation  region  shown  is  a  box 
which  contains  the  outer  section  of  a  rotor.  The 
boundary  conditions  on  the  faces  of  this  box  are 
usually  undisturbed  flow  conditions,  except  for 
the  inboard  face  where  two-dimensional  flow  is 
assumed.  The  surface  tangency  condition  is 
expressed  by  specifying  the  surface  slopes  on  a 
mean  plane  of  the  rotor.  In  order  to  represent 
vortices  and  vorticity  sheets,  interior  planes  of 
potential  discontinuity  are  Imposed.  In  Fig.  1, 
two  such  planes  are  shown.  One  plane,  which  is 
aft  of  the  trailing  edge,  models  the  trailing  and 
shed  vorticity.  Another  sheet  below  the  rotor 
models  a  vortex  located  at  the  edge  of  this  sheet. 
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The  model  shown  in  rig.  1  automatically  takes 
account  of  the  inflow  produced  by  the  vorticity  In 
the  computation  region.  However,  the  vorticity  in 
the  remainder  of  the  total  rotor  flow  field  is  not 
automatically  accounted  for.  The  simplest  way  to 
include  this  outer  flow  seems  to  be  to  compute  a 
partial  angle  of  attack,  a  ,  which  is  a  result  of 
the  vorticity  in  the  outer  region,  and  impose  it 
on  the  mean  surface  boundary  condition.  In 
return,  the  loads  computed  by  the  finite- 
difference  computation  provide  the  basic  vorticity 
distribution  for  the  cuter  flow  computation.  This 
exchange  of  flow  angle  and  load  information 
between  the  inner  and  outer  flow  regions  is  the 
heart  of  a  finite-difference  embedding  scheme. 

This  embedding  will  now  be  demonstrated  for  two 
different  flight  situations. 


Model  Hover  Testing 


The  simplest  problem  that  requires  the  above 
embedding  procedure  is  a  high-speed  hover. 

This  problem  is  relatively  simple  because  it  is 
steady,  and  the  required  validation  data  are 
fairly  easy  to  obtain  on  the  model  scale.  Fig¬ 
ure  2  shows  a  7.5-ft-diam  model  which  was  tested 
ir.  the  Aeromechanics  hover  facility.  The  model 
was  equipped  with  a  hub-mounted  Scanivalve  assem¬ 
bly  which  permitted  the  measurement  of  pressure  at 
'20  points  on  the  rotor.  Since  the  wake  geometry 
is  unpredictable,-  its  measurement  is  a  vital  ele¬ 
ment  of  such  testing.  These  measurements  were 
performed  with  traversing  hot  wires.  The  result¬ 
ing  information  was  then  used  as  the  required 
prescribed  wake  data  in  an  integral  flow  code — in 
this  case  the  Analytic  Method,  Inc.  code,  HOVER . 
For  subcrltical  tip  speeds,  HOVER  predicts  the 
measured  lift  distribution  very  well.  However, 
for  supercritical  tip  speeds,  HOVER  cannot  treat 
the  resulting  transonic  non! inearity ,  and  it 
becomes  necessary  to  use  an  embedded  finite- 
difference  solution.  For  the  present  problem  a 
steady  rotor,  small-perturbation  code  was  used.^ 
The  boundary  conditions  employed  were  the  same  as 
depicted  in  Fig.  ',  except  that  two  vortices  were 
included  in  the  computation  -egion.  In  the  course 
of  the  computation,  the  locations  of  tne  vortices 
were  fixed  by  the  experimental  data,  and  their 
strength  wa3  equated  to  the  maximum  blade  bound 
circulation.  A  typical  comparison  of  measured  and 
computed  pressure  distribution  is  3hown  in  Fig.  3 
for  a  tip  speed  W-r  -  0.377. 

The  fact  that  the  finite-difference  embedded 
solution  reproduces  the  measured  loads  well  beyond 
the  range  of  validity  of  the  linear  HOVER  code 
demonstrates  the  validity  of  the  small- 
perturbation  method.  It  also  shows  that  the 
scheme  keeps  consistent  account  of  vorticity,  in 


spite  of  this  vorticity  being  split 
entirely  separate  codes. 


Forward  Flight  Computations 


Transonic  flows  are  intrinsically  unsteady, 
and  advancing-olade  transonic  flows  can  be  among 
the  most  unsteady  flows  of  all.  Nevertheless, 
tabular  profile  data  take  no  account  of  this 
unsteadiness.  Such  flows  are  difficult  to  compute 
if  the  wake  is  to  be  accounted  for,  as  was  done  in 
the  previous  hover  case.  For  such  cases,  the 
vortices  from  previous  Slades  are  no  longer  neatly 
aligned  with  a  coordinate;  instead,  they  assume  a 
wide  range  of  orientations  as  tney  move  through 
the  grid.  At  present,  such  computations  are  being 
attempted  only  in  two  dimensions.  However,  a 
simple  first  approach  to  the  problem  is  to  account 
for  tne  inflow  entirely  with  a  partial  angle  of 
attack,  with  no  vortex  in  the  grid  other  than  the 
trailing  and  shed  vorticity  sheet.  Such  a  scheme 
is  almost  certainly  valid  at  high  advance  ratios 
where  the  induced  inflow  is  small.  This  scheme 
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has  been  implemented  in  a  code7  we  call  FDR 
(Finite-Difference  Rotor  Code);  FDR  solves  the 
low-frequency  variant  of  Eq.  (1)  (that  is,  A  *  0). 

Another  problem  with  forward  flight  computa¬ 
tions  is  that  the  blade  is  constantly  responding 
to  the  varying  loads,  and  a  trim  solution  must  be 
obtained.  The  problem  here  is  one  of  efficiency. 
Trim  solutions  are  iterative  processes  often 
requiring  the  rough  equivalent  of  computing 
'0  rotor  revolutions.  This  much  computing  is  no 
problem  when  the  blade  forces  are  obtained  from  a 
table,  but  it  would  be  prohibitive  for  a  finite- 
difference  solution.  A  solution  to  this  problem 
is  to  perform  the  basic  finite-difference  solution 
outside  the  central  trim  iteration.10  Inside  the 
trip  loop,  tables  (or  any  fast,  convenient  approx¬ 
imate  to  a  finite-difference  solution)  can  be  used 
to  find  a  load  correction  based  on  the  new  blade 
motion.  That  is. 


, ' ^old ' 


(a)  -  C 


u.  .  .  JiQ 

table 


where  a  and  are  the  angles  cf  attack  from 

-he  current  and  previous  trim  loops,  respectively. 
The  solution  is  converged  when  1  and  the 

lift  correction  vanishes;  t.nat  is,  when  the 
f inite-difTerence  lirt  is  filly  consistent  with 
the  rotor  .  of  lew  and  motion.  This  process  has 
been  followed  u3lng  the  comprehensive  -otor  code 
CAMRAD  to  provide  the  outer  inflow,  blade 
motion,  and  trim  computations. 

The  matching  of  the  1AMRAD  and  F‘P  c-'des  1  = 
summarized  in  the  •'low  iiagram  shewn  m  Fig.  « . 

The  process  is  started  by  obtaining  i  trimmed, 
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[  nonuniform  inflow  solution  with  the  lift  totally 

obtained  from  airfoil  tables.  This  Is,  In  fact, 
the  normal  operation  of  CAMRAD,  except  that  par¬ 
tial  wake-influence  coefficients  are  computed  In 
addition  to  the  usual  full-wake  values.  These 
partial  coefficients  are  used  to  obtain  partial 
,  angles  of  attack,  which  In  turn  are  fed  to  the 

I  FOR.  The  lift  distributions  so  obtained  are  then 

fed  to  the  CAMRAD  trim  loop  where  the  airfoil 
tables  are  used  to  find  a  lift  correction.  The 
process  now  iterates  between  the  FDR  and  the  trim 
loop  until  convergence  of  a  is  achieved.  This 
scheme  is  an  efficient  one  because  it  keeps  all 
the  time-consuming  computations  (the  influence- 
coefficient  and  especially  the  finite-difference 
calculations)  out  of  the  Innermost  trim  loop.  It 
has  been  our  experience  to  date  that  the  conver¬ 
gence  of  this  scheme  is  extremely  rapid. 

A  suitable  body  of  experimental  data  on  which 
to  test  the  present  code  was  obtained  at  a  recent 
model-rotor  test  performed  at  the  Deutsch- 
Nlederlaendischer  Windkanal  (DNW).1^  The  rotor  Is 
a  lf7-scale  model  of  the  Cobra  operational  load 
survey  (OLS)  blades  (Fig.  5).  -ike  the  OLS 
blades,  the  model  Is  pressure  Instrumented  (but 
not  nearly  so  densely).  At  the  present  time,  a 
small  amount  of  the  pressure  data  is  available. 
Moreover,  the  model  Is  fairly  uncomplicated  and  la 
a  good  first  test  case  for  the  present  method. 

The  rotor  Is  a  teetering  rectangular  blade  with  a 
linear  twi3t  and  has  a  modified  BHC-5<I0  airfoil. 
Since  the  model  Is  quite  stiff,  we  have  not  con¬ 
sidered  blade  bending  or  twisting  in  the  present 
computations. 

The  available  data  sets  are  high-speed  cases. 
A  typical  run,  which  we  have  computed,  has  a  tip 
Mach  number  Mj  of  0.663,  an  advance  ratio  u  of 
0.298,  and  a  C^/o  •  0.0769.  For  this  case,  the 
contribution  of  the  grid-enclosed  trailing  vortex 
to  the  full  angle  of  attack  i3  seen  in  the  span- 
wise  distributions  of  a  and  a  at  *  *  90° 

(Fig.  6).  At  this  azimuth,  a  ?ip  vortex  passes 
beneath  the  rotor  at  0.7R.  Accordingly,  a  is 
augmented  outboard  and  diminished  inboard  of  this 
location.  Since  the  rotor  lift  varies  similarly, 
the  trailing  vorticity  will  have  an  opposing 
effect  on  the  Induced  flow  of  the  vortex.  Remov¬ 
ing  this  opposing  effect  results  in  an  a  dls- 
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trlbution  with  a  much  stronger  variation.  This 
large  difference  between  a  and  a  demonstrates 
the  importance  of  keeping  a  consistent  account  of 
the  wake  with  a  matched  scheme.  In  Fig.  6,  the 

distribution  of  a,  a  ,  C,  (lift  computed  in 
P  LFD 

FDR),  and  C^  (the  lift  computed  in  the  CAMRAD 
trim  loop)  are  shown  for  three  iterations.  The 
N  ■  9  iteration  is  the  starting  CAMRAD  lifting- 
line  solution.  All  subsequent  iterations  are 
complete  Integral -differential  cycles.  For  N  *  ' 


there  is  a  moderate  difference  between  the  C, 

“FD 

and  the  C^.  However,  for  N  *  2  there  is  no 
difference  between  them,  and  the  solution  appears 
completely  converged.  Another  view  of  the  conver¬ 
gence  is  shown  in  Fig.  7,  which  shows  the  azi¬ 
muthal  variation  of  a  for  various  iterations. 
Again  the  solution  appears  completely  converged  by 
the  second  iteration,  with  the  greatest  change 
occurring  In  the  second  quadrant. 

Figure  8  shows  the  spanwlse  variations  of 
C^  for  a  wide  range  of  azimuths  for  the  pure 
CAMRAD  and  the  converged  CAMRAD-FDR  solutions.  It 
is  seen  here  that  the  effect  of  switching  on  the 
finite-difference  solution  is  to  shift  the  loading 
somewhat  from  the  rear  to  the  front  of  the  rotor 
disk.  Accordingly,  the  greatest  trim  erfect  Is  on 
the  lateral  cyclic  input  (Fig.  9).  Furthermore, 
the  effect  of  the  finite-difference  solution  is  to 
bring  the  trim  settings  into  somewhat  closer 
accord  with  the  experimental  values.  Note  also 
that  the  trim  settings  are  essentially  fixed  after 
the  first  Iteration. 

At  present,  the  only  available  pressure  data 
are  for  the  upper  blade  surface  at  r/R  -  0.95. 

The  comparison  of  the  measured  and  computed  uppei — 
surface  pressures  is  shown  In  Fig.  10.  A  similar 
pressure  comparison  at  an  advance  ratio  of  0-3*15 
is  shown  In  Fig.  11.  There  Is  a  considerable 
difference  between  measured  and  computed  pressures 
In  the  first  quadrant  where  the  pressures  appear 
to  be  generally  underpredicted.  The  comparison 
tends  to  be  better  in  the  second  quadrant,  where 
the  shock  location  is  correctly  predicted. 


Concluding  Remarks 


The  above  results  are  an  indication  of  the 
most  probable  immediate  direction  of  finite- 
difference  computations  for  application  to  high¬ 
speed  rotor  loads.  The  mixed  integral- 
differential  approach  to  transonic  load  prediction 
has  the  advantage  of  making  immediate  use  of  our 
existing  loads  codes.  It  also  can  greatly  reduce 
the  need  for  two-dimensional  airfoil  testing.  (It 
should  be  borne  In  mind  that  two-dimensional  pro¬ 
file  data  are  probably  of  dubious  value  when  non- 
rectangular  planforms  or  transonic  unsteadiness 
are  involved.)  In  its  present  stage  of  develop¬ 
ment,  the  approach  gives  fair  to  good  comparison 
with  experimental  pressure  data.  However,  consid¬ 
erably  more  comparisons  with  forward  flight  data 
are  required.  The  method  clearly  can  be  enhanced 
by  employing  a  conservative,  full-potential 
f ini te-dlf ference  model  and  by  improving  the  neai — 
wake  modeling  capability.  This  latter  improvement 
would  involve  inserting  the  moving  vortex  back 
into  the  finite-difference  grid.  Such 
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modifications  will  certainly  be  required  if  the 
method  is  to  be  extended  to  lower  advance  ratios 
at  which  the  blade/vortex  interactions  are  more 
important .  Even  in  its  present  form,  however,  the 
mixed  integral-differential  scheme  affords  a 
considerable  Improvement  in  our  flow  prediction 
capabi 1 i t ies . 

In  the  future,  finite-difference  modeling  can 
be  developed  and  extended  for  application  to  a 
wider  variety  of  problems,  the  most  important  of 
wntcn  is  probably  the  prediction  of  rotor-wake 
behavior.  Ultimately,  it  will  not  be  unreasonable 
to  expect  the  development  of  unified  finite- 
difference  models  of  the  entire  rotor/fuselage. 
However,  this  goal  will  require  many  fundamental 
advances  in  numerical  modeling  and  in  basic  flow 
physics. 
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